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Motivated by the observation of unusual magnetism in Cea;Cu2Si2 (a;~ 1), we study the effect of 
disorder, such as Ce vacancy, on the competition between superconductivity (SC) and antiferro- 
magnetism (AF) on the basis of the phenomenological Ginzburg-Landau theory. Assuming that the 
AF-SC transition is of first order in clean system, we show that a single impurity in the SC state 
can induce staggered magnetization by suppressing the SC around it. For finite concentration of 
impurities, the first-order AF-SC boundary in the clean case is replaced by a finite region where the 
SC and the induced AF moments coexist microscopically with spatially varying order parameters. 
We argue that spin excitation spectrum in the coexistent state has a dual structure of SC gapped 
mode and the low-energy spin-wave mode. In accordance with the emergence of AF out of SC 
ground state, the spectral weight will be transferred from the former mode to the latter, keeping 
the structure of both modes basically unchanged. 



I. INTRODUCTION 

The interplay between magnetism and superconductiv- 
ity (SC) is one of the interesting and profound phenom- 
ena ubiquitous in strongly correlated electron systems 
such as high-T c cuprates, 'low- TV ruthenates, heavy- 
electron compounds and molecular solids. Especially, 
heavy-electron systems often exhibit mysterious mag- 
netism and exotic superconductivity, most of which re- 
main to be characterized |jj . 

The heavy-electron compound Cej;Cu2Si2 (x ~ 1) is 
located near the phase boundary between superconduc- 
tivity and antiferromagnetism (AF), and shows an ex- 
treme sensitivity to the variation of the Ce concentration 
x l^-fl- Upon reducing x, the ground state is trans- 
formed from the spin-singlet SC (x = 1.025, 1.00) to AF 
(x = 0.975). For x = 0.99, several experiments indicate 
the coexistence of SC and the so-called A-phase, the lat- 
ter being characterized by the presence of spin fluctua- 
tions with unusually low frequency. Recent NQR mea- 
surements by Ishida et al. M revealed that the energy 
scale of spin fluctuations in the A-phase is extremely low 
and of the order of NQR frequency, Wnqr ~1-10 MHz. 
Zero-field /.(SR studies also indicated the existence of ran- 
dom or slowly- fluctuating internal magnetic field ||-p|. 

In this letter, we propose that in Cea;Cu2Si2 with such 
a particular type of magnetism an essential role is played 
by disorder JljJ. A Ce vacancy, for example, can be the 
origin of such disorder, which will strongly disturb the 
electron coherence around it. In the SC state, it will act 
as a pair breaker and locally reduce the SC amplitude. It 
may even induce local magnetization, possibly the stag- 
gered one, around it. We present a phenomenological de- 
scription based on the Ginzburg-Landau (GL) theory for 
such effects of disorder on the SC state. A key assump- 
tion made is that the SC state is in competition with 



the AF phase with a first-order phase boundary in clean 
case, which enables the AF phase to nucleate where SC 
is suppressed. We limit the discussion to absolute zero, 
namely, the ground state. 

Disorder-induced AF in spin-gap systems has been rec- 
ognized through studies of impurity effects in inorganic 
spin-Peierls system CuGeOa [ll],[L2[. In this system, a 
small amount of replacement of Ge by Si, or Cu by Zn 
induces local staggered moment with local reduction of 
lattice dimerization, leading to the coexistence of both 
long-range orderings (LRO's), lattice dimerization and 
AF, since there is no frustration for AF as emphasized in 
ref . p2] . The present work attempts a general treatment 
for such phenomena, especially in higher dimensions, us- 
ing the GL theory. Details will appear elsewhere. 



II. CLEAN SYSTEMS 

We begin by considering a clean system, i.e., without 
impurities, described by the GL free energy density 
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Both order parameters, i/j and M, for SC and AF, respec- 
tively, are assumed to be real scalars |1J. We are con- 
cerned only with their amplitude modulations here and 
neglect the phase (orientational) degrees of freedom for 
ip (M). The constants, a, (3, ccm, Pmi 7j P an d Pm, are all 
assumed to be positive. Especially, we consider the case 
where the ground state is always in a broken-symmetry 
phase of either SC or AF, which is expressed by the condi- 
tions, a > and cum > 0. This viewpoint has been taken 
by Zhang in the SO(5) theory Ej]. We also assume that 
there is no coexistence of SC and AF in the clean case, 



which is ensured by the condition, 7' > minl/3',/3^}, 
where 7' = j/aatM, P' = fl/oi 2 and /3' M = Pm/c^m- Then 
the ground state is SC if j3' < (3' M and AF otherwise. 

In the following, we consider the parameter region 
where the SC state is stabilized in the bulk clean sys- 
tem, but is located near the phase boundary to the AF 
state. The SC order parameter takes a constant value 
tpo = y/a/(3 in the bulk in this case. 



III. SINGLE IMPURITY 

We now introduce a single impurity (at the origin) into 
the otherwise uniform SC ground state. For CeCu2Si2, 
this may represent a Ce vacancy, for example. Assuming 
that the impurity acts as a local pair breaker for SC and 
locally favors AF relative to SC, we model its effect by 



a — ► a(r) =a[l-V 6(r)], 
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where Vo, Vi > 0. The GL equations are given by 
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Equation (Q) describes the suppression of the SC order 
parameter ip around the impurity. If we neglect the 7 
term for simplicity [lq| , the solution is written as 



*!>(?) = Vojpf(r ) Z,Vo) 



(6) 



with £ = yjp/2a and r = |r|. Qualitative features of 
the function / are as follows. It is suppressed apprecia- 
bly from the bulk value of 1 only within a region r <^£, 
where it decreases almost linearly with r as r — ► 0. The 
value f(r = 0) at the origin is a monotonically decreasing 
function of Vo; it is 1 for Vq — and vanishes as Vo — > 00. 
Equation (||) describes the possible emergence of stag- 
gered moment M around the impurity. To study this, we 
first drop the non-linear term (~ Af 3 ) and consider the 
following eigenvalue problem 



-^V 2 + C/ M (r) 



M (r) = eM(r), 



where the 'potential' Um (r) is given by 
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In eq.(ra), the first term, 6cm = ~ a M + tV'oj i s positive 
p6[, which assures the local stability ('local' in the if>-M 
space) of the pure SC state in the clean case. The sec- 
ond and the third terms represent 'attractive potential' 
due to the suppression of ip, and the direct coupling to 
the impurity, respectively. If we have a non-zero solution 



M(r) with negative 'energy' e, this signals the appear- 
ance or 'nucleation' of the AF phase around the impu- 
rity in the otherwise SC system ]r7p . Its spatial pattern 
is given by the corresponding eigenfunction (with lowest 
eigenvalue, e) and the amplitude is determined by bal- 
ancing the quadratic term (~eM 2 ) with the quartic term 
(~/3m M 4 ), which leads to M ~ y/-e//3 M . 

As we have seen, the characteristic length-scale for ip 
is given by £. On the other hand, spatial variation of M 
occurs on the scale, £m = \fPMj2otM. Q Depending on 
their ratio, £/£m, there are two different regimes, namely, 
£ "C S,m (type I) and £ 3> £m (type II). The intermedi- 
ate regime (£ ~ £ M ) may be the case where the SO (5) 
symmetry, a rotational symmetry in the tp-M space, is 
expected. We examine these three cases in the following. 

(i) Type I (£ <C £m)'- The suppression of tp is limited 
to a small region compared to £m- Its effect on M can 
be shown to be negligible, and the Vi-term is necessary 
for M to nucleate. Rewriting eq.(R) as 



(-V 2 + K 2 M )M(r) = VxModiv), 



(9) 



with v\ = Vi/£ 2 M , km = «m(e) = V 2 ((5m - £)/pm, and 
Mq = M(r = 0), the solution is obtained as 



M(r)= Vl M G{r;K M ), 



(10) 



where G(r; k) is the Green's function satisfying (—V 2 
K 2 )G(r; k) — 5(r). The explicit form is given by 



G(t;k) 
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r /4irr, (11) 



for dimensions d = 1,2 and 3, respectively, where _ftTo is 
the modified Bessel function of rank zero. (The present 
Ce-compound corresponds to d = 3.) Equation (10|) ex- 
presses a nucleated staggered magnetization if the solu- 
tion exists such that Mq ^ and e < 0. 

The condition for Mq ^ is obtained from ( |l0[ ) by 
putting r = and dividing both sides by M as 



1 = v 1 G{r = a\n M {e)). 



(12) 



This is an eigenvalue equation for e. For d = 2 and 3, 
the right hand side is evaluated by introducing a cut-off 
a which will be of the order of the lattice constant [p| . 
The nucleation of M occurs for v\ > V\ tCT where 

ui.cr = 2/|m, 2ir/K (a/^M), 47raexp(a/| A f), (13) 



for d = 1, 2 and 3, respectively, with £m = y/ pm/^olm- 

(ii) Type II (£ ;» £m): In this case, -0 is suppressed in 
a wider region compared to £m, and it is possible for M 
to nucleate even without the Vi-term. As an example we 
consider the case Vo — > 00 and V\ = 0. Then the potential 
for M is given by a harmonic one generated by the linear 
r-dependence of ip( r ) around the impurity. If we assume 
such harmonic potential in the entire space, the condition 
for the nucleation of M is obtained as £m < £/d 2 . 



(iii) '50(5) symmetric' case (£ = £m): In case 
both length-scales are comparable, we must solve the 
coupled equations for i/> and M simultaneously. How- 
ever, a simplification is possible if the 'anisotropy' is 
weak in the ip-M space. We define the case what may 
be called '50(5) symmetric' by £ = £m- Wc then 
introduce the angular (9) and radial (R) variables by 
(a 1//2 V>, Um M) = (Rcos9,Rsm9), to obtain 
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for clean systems, where [5(6) = [5 + g cos 29 + B sin 2 29, 
00 = W+ Pm)> 9 = W ~ P'm), and B = ±( 7 ' - ft). 
The terms with g and B represent, respectively, the 
'anisotropy energy' and the energy barrier between SC 
and AF. (g < favors SC relative to AF.) The con- 
dition of weak anisotropy, \g\, \B\ <C [5q, restricts the 
low-energy configurations to the manifold R = R(9) = 
[5(9)^ 1 / 2 , which as a function of 9 has small variation 

— 1/2 

compared to the average ~ Rq = p . Introducing 
the impurity potential by F; mp = VR 2 cos 296 (r), and 
considering the spatial variation of the order parame- 
ters, we again recognize two length-scales, £r ~ £ and 
£g ~ £y/[5o/max{\g\, \B\}, now associated with variables 
R and 9, respectively. From the above condition, we have 
£r <C £g. Therefore, while both R and 9 deviates from 
the bulk values in the vicinity of the impurity (r < £r), 
for r > £r, only 9 governs the spatial variation and R 
follows 9 completely as R = R(9) given above. The be- 
haviours in the region r > £r will hence be described by 
the effective free energy, Fg: 

Fg/Rl =£ 2 (V9) 2 +g(r)cos29 + B'sm 2 29, (15) 

where g(v) = (g/2/5 ) + V6(r), B' = B/2[5„. For B = 
and d = 1, the present problem is equivalent to that of 
the disordered spin-Peierls system |1J,|2(J ■ The one with 
B = and d = 2 (with V = 0) corresponds to the 50(5) 
description proposed for high-T c cuprates Jl4[ . 



IV. FINITE CONCENTRATION OF IMPURITIES 

In the presence of many impurities (located at {Ri}), 
the solution is written as 



M(r)=^m i M 1 (r-R i ) 



(16) 



in the dilute concentration limit. Here M\ is the solution 
for the single-impurity problem obtained in the previous 
section, and we have recovered the vector nature of the 
AF order parameter. The unit vector, nii, represents the 
direction of the staggered moment nucleated around the 
impurity at Rj, and the interaction among them will be 
described by 



i<j 



The effective exchange constants 
Jij = (-2e) /drAfi(r - R,)Mi(r 
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(17) 



(18) 



do not have frustration in the present model focussing 
on the commensurate AF, but will allow it in general. 
For depleted quantum spin systems, similar models have 
been obtained previously fl2lf . 

Orientational fluctuations of mi's will produce low- 
energy spin excitations even when the LRO of SC is not 
completely destroyed by the emergence of M. If mj's 
align and exhibit AFLRO, their fluctuations will show up 
as a gapless Goldstone (spin-wave) mode. We note that 
the coexistence of both types of LRO, SC and impurity- 
induced AF, is possible due to the spatial variation of 
the two order parameters in spite of their competition. 
This general aspect has been first noticed in studies on 
disordered spin-Peierls systems. |13] If AFLRO failed to 
be established for some reasons (see below), fluctuations 
of mi's will remain as a diffusive mode with very low en- 
ergy with the total spectral weight being proportional to 
M 2 and the concentration of impurities. 



V. DISCUSSION 

We have given an analysis based on the GL theory of 
the effect of disorder in the SC ground state which is in 
competition with the AF state. Our interest is focussed 
on the nucleation of the AF phase around disorder cen- 
ters. As the concentration of impurities is increased, the 
original SC state will first change into the microscopically 
coexistent state with nucleated staggered moments, and 
finally be transformed into the complete AF. The first- 
order AF-SC boundary in the clean case thus turns into 
a finite region of coexistence. 

Analogously to disordered spin-Peierls systems |g2|,^3[ , 
the present coexistent state has spin excitations with 
two distinct features, the SC gapped excitation due to 
quasiparticles and the gapless spin-wave mode, well sep- 
arated from each other by a 'transparent' region (see 
Fig.pl). In accordance with the emergence of AF out of 
SC ground state, the spin excitation spectrum will change 
mainly through the transfer of the spectral weight from 
the gapped mode to the spin-wave mode, keeping the 
structure and energy scales of both modes essentially un- 
changed, with a possible reduction of the SC gap and an 
increase of the spin-wave velocity. 
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FIG. 1. Spin excitation spectrum in the coexistent phase 
of SC and disorder-induced staggered magnetization. It has a 
dual structure of SC gapped mode from quasiparticles and the 
low-energy mode from nucleated staggered moments, (a) In 
the absence of AFLRO, low-energy excitations are diffusive, 
(b) In the presence of AFLRO, low-energy excitations are 
Goldstone (spin wave) modes. The dot at Qaf indicates the 
Bragg spot of AF. The energy scale for the low-energy mode 
is exaggerated. 

We would like to propose that the above picture de- 
scribes well the evolution of the ground state from SC 
to AF in Ce x Cu2Si2. Especially, we propose that the A- 
phase realized near the phase boundary is characterized 
by a coexistence phase of SC and the nucleated staggered 
moments. Experiments so far performed on samples in 
the A-phase did not detect any signals of magnetic LRO 
p4J , but revealed the existence of dynamical fluctuations 
with very low frequency of NQR scale ||. Such low- 
energy fluctuations will be naturally understood as orien- 
tational fluctuations of the nucleated staggered moments 
around impurities. In the absence of true AF LRO, they 
behave as diffusive (Fig. 1(a)). The spectral weight can 
be substantial at low energy since they are transformed 
into a Bragg spot and collective spin-wave modes as AF 
LRO is reached (Fig.l(b)). 

The nucleation of staggered moments will occur at 
some finite temperature. Then what is their fate in the 
low-temperature limit where thermal fluctuations are di- 
minished? There are several possibilities depending on 
the interaction Jy among them; (i) Constructive interac- 
tion (without frustration) leads to AFLRO, (ii) the pres- 
ence of frustration leads to a spin-glass state, (iii) even 
if there is no frustration, the interaction among them is 
so weak and their orientations freeze almost randomly 
and independently, resulting in a spin-glass-like state, 
and (iv) they continue fluctuating due to quantum ef- 
fects. The spin-Peierls system Cui_ x Zn x Ge03 offers an 
example of case (i), where AF LRO has been observed in 
samples with x being as low as 1.12 x 10 -3 with Neel 
temperature T N = 0.0285K. Q For Ce x Cu 2 Si 2 with 
x = 0.99, the NQR experiment has been done down to 
temperature T = 0.012K, which suggests the case (ii) or 
(iii). If the nucleated staggered moments are incommen- 
surate [ pJ6| , their interaction will be frustrated, leading 
to the case (ii). 

Our present picture predicts the spectral weight trans- 
fer between the 'rigid bands' of gapped mode and the 



'spin-wave' mode (Goldstone or diffusive, depending on 
the presence or absence of AF LRO) as the ground state 
is tuned from SC to AF through the coexistent phase. It 
will be interesting to explore experimentally this spectral- 
weight transfer by sweeping through the phase transition 
via, for example, pressure as proposed in ref. M. 

In the coexistent state proposed in this paper, SC and 
AF are in competition even though they 'coexist'. The 
/LtSR studies ||[JJ seem to have revealed their compet- 
ing nature experimentally, but then concluded a spatial 
seggregation of the two phases. Before deriving such a 
conclusion, an analysis from the present viewpoint may 
be necessary. 

In conclusion, we have examined a possible nucleation 
of the AF phase due to disorder in the otherwise pure 
SC state and proposed that such a coexistent phase is 
responsible for the peculiar magnetism in the so-called 
A-phase of Ce^Cu2Si2. Such effects of disorder on the 
interplay between magnetic ordering and superconduc- 
tivity (or spin gap in general) will be quite ubiquitous 
in strongly correlated systems and deserve attention in 
interpreting experimental data. 

We would like to thank Professor Y. Kitaoka for a nice 
introduction to Ce x Cu2Si2, as well as for valuable com- 
ments. 
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